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ABSTRACT 



^ ■ The unusual supernova SN1998bw, which is thought to be related to the 7- 

■ ray burster GRB980425, is a possible link between the two classes of objects. 

Analyzing the extensive radio emission data available for SN1998bw, we are able 
to describe its time evolution within the well established framework available for 
the analysis of radio emission from supernovae. This then allows description of 
a number of physical properties of the object. The radio emission can best be 
explained as interaction of a mildly relativistic (F ~ 1.6) shock with a dense pre- 
explosion stellar wind established circumstellar medium (GSM) which is highly 
structured both azimuthally, in clumps or filaments, and radially, with two ob- 
served density enhancements separated by ~ 3 x 10^^ cm. With assumptions 
as to pre-explosion stellar wind conditions, it is possible to estimate that the 
progenitor to SN 1998bw had a mass loss rate of ~ 3.5 x 10"^ M0 yr~^ with 



Also Astrophysics Division, Space Science Department of ESA. 



- 2 - 



at least two ~ 30% increases in mass-loss rate; the most recent extending from 
~ 1, 600 — 4, 700 yr before explosion and the oldest known having occurred, with 
possibly comparable length, ~ 12, 000 yr before explosion. Because of its unusual 
characteristics for a Type Ib/c SN, the relation of SN1998bw to GRB980425 is 
strengthened with consequent improvement in our understanding of these poorly 
understood objects. 



1. Introduction 

While generally accepted that "most" GRBs are extremely distant and energetic (see, 
e.g., Paczynski 1986; Goodman 1986), the discovery of GRB980425 (Soffitta et al. 1998) 
on 1998 April 25.90915 and its possible association with a bright supernova, SN1998bw 
at RA(J2000) = 19''35'"03!31, Dec(J2000) = -52°50'44'.'7 (Tinney et al. 1998), in the 
relatively nearby spiral galaxy ESO 184-G82 aX z = 0.0085 (distance ~ 40 Mpc for Hq = 65 
km s-^ Mpc^i) (Galama et al. 1998a,b; Lidman et al. 1998; Tinney et al. 1998; Sadler 
et al. 1998; Galama a/. 1999; Woosley, Eastman, & Schmidt 1999), has introduced the 
possibihty of a SN origin for at least some types of GRBs. The estimated explosion date of 
SN1998bw in the interval 1998 April 21 - 27 (Sadler et al. 1998) corresponds rather well 
with the time of the GRB980425 outburst. Iwamoto et al. (1998) feel that they can restrict 
the core collapse date for SN1998bw even more from hydrodynamical modeling of exploding 
C + O stars and, assuming that the SN1998bw optical light curve is energized by ^^Ni decay 
as in Type la SNe, they then place the coincidence between the core collapse of SN1998bw 
to within +0.7/-2 days of the outburst detection of GRB980425. 

Classified initially as an SN optical Type lb (Sadler et al. 1998), then Type Ic (Patat & 
Piemonte 1998), then peculiar Type Ic (Kay, Halpern, & Leighly 1998; Filippenko 1998), 
then later, at an age of 300 - 400 days, again as a Type lb (Patat et al. 1999), SN1998bw 
presents a number of optical spectral peculiarities which strengthen the suspicion that it 
may be the counterpart of the 7-ray burst. 

When the more precise BcppoSAX NFl was pointed at the BcppoSAX error box 10 
hours after the detection of GRB980425, two X-ray sources were present (Plan et al. 1999). 
One of these, named SI by Pian et al. (1999), is coincident with the position of SN 1998bw 
and declined slowly between 1998 April and 1998 November. The second X-ray source, S2, 
which was ~ 4' from the position of SN 1998bw, was not (or at best only marginally with 
a < 3(7 possible detection six days after the initial detection) detectable again in follow up 
observations in 1998 April, May, and November (Pian et al. 1999, 2000). 
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However, although a bit of concern remains that the Pian et al. (1999, 2000) X-ray 
source, S2, might have been the brief afterglow from GRB980425 rather than the Pian source 
SI associated with SN1998bw, Pian et al. (2000) conclude that their source SI has a "high 
probability" of being associated with GRB980425 and that S2 is more likely a variable field 
source. 

In order to analyze the radio emission from SN1998bw we shall first collect, rationalize, 
and extend the formalism which has been so successful in understanding the radio emission 
from radio supernovae (RSNe) (Section 2), and apply it to SN1998bw (Section 3). From 
this we are able to obtain a detailed description of the radio light curves for the object 
and physical parameters of the SN system (Section 4). Finally, we compare a number of 
properties of Type Ib/c SNe and GRBs (Section 5). 

2. Radio Supernovae 
2.1. Background 

A series of papers pubhshed over the past 20 years on radio supernovae (RSNe) has 
established the radio detection and, in a number of cases, radio evolution for approximately 
two dozen objects: 3 Type lb supernovae (SNe), 3 Type Ic SNe (N.B.: Because the differences 
between the SN optical classes are slight - Type lb show strong He I absorption while Type 
Ic show weak He I absorption - and there are no obvious radio differences, we shall often 
refer to the classes as Type Ib/c), and the rest Type II SNe. A much larger list of more 
than 100 additional SNe have low radio upper limits and can be found at: 
http://rsd-www.nrl.navy.mil/72i4/weiler/kwdata/RSNtable.txt. 

In this extensive study of the radio emission from SNe, several effects have been noted: 
1) No Type la SN has ever been detected at radio wavelengths; 2) Type Ib/c SNe are radio 
luminous with steep spectral indices (generally a < —1; S oc i>~^") and a fast turn-on/turn- 
off, usually peaking at 6 cm near or before optical maximum; and 3) Type II SNe show a 
range of radio luminosities with flatter spectral indices (generally a > — 1) and a relatively 
slow turn-on/turn-off, usually peaking at 6 cm significantly after optical maximum. Type 
Ib/c may be fairly homogeneous in some of their radio properties while Type II, as in the 
optical, are quite diverse. 

Measurements of the RSN multi-frequency radio light curves and their evolution with 
time show the density and structure of the CSM, evidence for possible binary companions, 
dumpiness or filamentation in the presupernova wind, mass-loss rates and changes therein 
for the presupernova stellar system and, through stellar evolution models, estimates of the 
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ZAMS presupernova stellar mass and the stages through which the star passed on its way 
to explosion. A summary of the radio information on SNe can be found at: 
http:/ / rsd-www. nrl. navy, mil/ 7214 / weiler/ sne-home. html. 

2.2. RSN Models 

All known RSNe appear to share common properties of: 1) nonthermal synchrotron 
emission with high brightness temperature; 2) a decrease in absorption with time, resulting 

in a smooth turn-on first at shorter wavelengths and later at longer wavelengths; 3) a power- 
law decline of the flux density with time at each wavelength after maximum flux density 
(optical depth ~ 1) is reached at that wavelength; and 4) a flnal, asymptotic approach of 
spectral index a {S (X z/+") to an optically thin, nonthermal, constant negative value (Weiler 
et al. 1986; Weiler, Panagia, & Sramek 1990). 

Chevalier (1982a,b) has proposed that the relativistic electrons and enhanced magnetic 
fleld necessary for synchrotron emission arise from the SN blastwave interacting with a 
relatively high density CSM which has been ionized and heated by the initial UV/X-ray flash. 
This CSM is presumed to have been established by a constant mass-loss (M) rate, constant 
velocity (wwind) wind {i.e., p oc ^ ^2 ) from a massive stellar progenitor or companion. 
This ionized CSM is the source of some or all of the initial absorption although (Chevaher 
1998) has proposed that synchrotron self- absorption (SSA) may play a role in some objects. 

A rapid rise in the observed radio flux density results from a decrease in these absorption 
processes as the radio emitting region expands and the absorption processes, either internal 
or along the hne-of-sight, decrease. Weiler, Panagia, & Sramek (1990) have suggested that 
this CSM can be "clumpy" or "filamentary," leading to a slower radio turn-on, and Montes, 
Weiler, & Panagia (1997) have proposed at least one example for the presence of a distant 
ionized medium along the line-of-sight which is time independent and can cause a spectral 
turn-over at low radio frequencies. In addition to clumps or filaments, the CSM may be 
structured with significant density irregularities such as rings, disks, shells, or gradients. 

2.3. RSN Pcirameterized Radio Light Curves 

Following Weiler et al. (1986), Weiler, Panagia, & Sramek (1990), and Montes, Weiler, 
& Panagia (1997), we adopt a parameterized model (N.B.: The notation is extended and 
rationalized here from previous publications. However, the "old" notation of r, r', and r", 
which has been used previously, is noted, where appropriate, for continuity.): 
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5(mJy) = K, ( — ^) ( ^ ) e--" ( ^ ) ( ^ ) (1) 

2.3.1. Eocternal Absorption: Uniform & Distant 

^external = ''"CSMuniform + ''"distant = T + t" (2) 



rcsM.„,_ - r - j [j^] (3) 



ra.tant = r" = X4(^)"' (4) 



^.5.^. External Absorption: Clumpy or Filamentary 



TCSM._ =t' = K, (^) (5) 

with Xi, ^^2, -f^3, and corresponding, formally, to the flux density (i^i), uniform (^2, 
Ki), and clumpy or filamentary {K^) absorption, at 5 GHz one day after the explosion 
date to- The terms TcsMuniform ''"csMdumps describe the attenuation of local, uniform CSM 
and clumpy CSM that are near enough to the SN progenitor that they are altered by the 
rapidly expanding SN blastwavc. TcsMuniform produced by an ionized medium that uniformly 
covers the emitting source ("uniform external absorption"), and the (1— e~^^^^*<:i™p")rQg^j^j ^ 
term describes the attenuation produced by an inhomogeneous medium with optical depths 
distributed, with equal probability, between and TcsMdumps ( "clumpy absorption" ; Sec Natta 
& Panagia 1984 for a more detailed discussion of attenuation in inhomogeneous media). The 
^distant term describes the attenuation produced by a homogeneous medium which uniformly 
covers the source but is so far from the SN progenitor that it is not affected by the expanding 
SN blastwave and is constant in time. All absorbing media are assumed to be purely thermal, 
singly ionized gas which absorbs via frcc-frcc (f-f) transitions with frequency dependence z/~^'^ 
in the radio. The parameters 8 and 5' describe the time dependence of the optical depths 
for the local uniform and clumpy or filamentary media, respectively. 
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The f-f optical depth outside the emitting region is proportional to the integral of the 
square of the CSM density over the radius. Since in the simple Chevalier model the CSM 
density decreases as r~^, the external optical depth will be proportional to r~^, and since 
the radius increases as a power of time, r oc with m < 1 {i.e., m = 1 for undecelerated 
blastwave expansion), it follows that the deceleration parameter, m, is 

m = -S/3. (6) 

The model by Chevalier (1982a,b) relates /? and S to the energy spectrum of the relativistic 
particles 7 (7 = 2q; — 1) by 5 = a — /3 — 3 so that, for cases where K2 — and S is, therefore, 
indeterminate, we can use 

m = -(a-/?-3)/3. (7) 



2.3.3. Internal Absorption: SSA & Mixed f-f Absorption/Nonthermal Emission 

Since it is physically realistic and may be needed in some RSNe where radio observations 
have been obtained at early times and high frequencies, we have also included in Eq. 1 the 
possibility for an internal absorption term^. This internal absorption (rinternai) term may 
consist of two parts - synchrotron self- absorption (SSA; TintemaissA)' mixed, thermal f-f 
absorption/non-thermal emission (rinternaiff)- 

''"internal ''"internalssA ~^ ^intcrnalff (8) 



S" 



TinternalssA " ^^5 GHz J V 1 day J 



with and Kq corresponding, formally, to the internal, non-thermal (i/""^-^) SSA (K^) and 
the internal, mixed with nonthermal emission, thermal {j^^^'^) f-f absorption (Kq), respec- 



^Note that for simplicity we use an internal absorber attenuation of the form ( — "'^"■ntemai \ ^ 
is appropriate for a plane-parallel geometry, instead of the more complicated expression (e.g., Osterbrock 
1974) valid for the spherical case. The assumption does not affect the quality of our analysis because, to 
within 5% accuracy, the optical depth obtained with the spherical case formula is simply three-fourths of 
that obtained with the plane-parallel slab formula. 
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tively, at 5 GHz one day after the explosion date to- The parameters 5" and 5"' describe the 
time dependence of the optical depths for the SSA and f-f internal absorption components, 
respectively. 

A cartoon of the expected structure of an SN and its surrounding media is presented in 
Fig. 1 (see also, Lozinskaya 1992). The radio emission is expected to arise near the blastwave 
(Chevaher & Pransson 1994). 

2.4. RSN Results 

The success of the basic parameterization and model has been shown in the relatively 
good correspondence between the model fits and the data for all subtypes of RSNe: e.g., 
Type lb SN1983N (Sramek, Panagia, & Weiler 1984), Type Ic SN1990B (Van Dyk et al. 
1993), Type II SN1979C (Weiler et al. 1991, 1992a; Montes et al. 2000) and SN1980K 
(Weiler et al. 1992b; Montes et al. 1998), and Type Iln SN1988Z (Van Dyk et al. 1993; 
Lacey et al. 2001). (Note that after day ~ 4000, the evolution of the radio emission from 
both SN1979C and SN1980K deviates from the expected model evolution, that SN1979C 
shows a sinusoidal modulation in its flux density prior to day ~4000, and that SN1988Z 
shows a double drop in its flux density, after day ~ 2400 and again after day ~3100.) A 
more detailed discussion of radio observations of SNe and of modeling results is given in 
Weiler et al. (2001) 

Thus, the radio emission from SNe appears to be relatively well understood in terms 
of blastwave interaction with a structured CSM as described by the Chevaher (1982a,b) 
model and its modiflcations by Weiler et al. (1986), Weiler, Panagia, & Sramek (1990), 
and Montes, Weiler, & Panagia (1997). For instance, the fact that the uniform external 
absorption exponent 5 is ~ —3, or somewhat less, for most RSNe is evidence that the 
absorbing medium is generally an p oc r^^ wind as expected from a massive stellar progenitor 
which explodes in the red supcrgiant (RSG) phase. Also, the relatively smooth behavior of 
the radio light curves suggests that both the external medium and the relativistic electron 
acceleration process are rather regular and smooth. In addition, there are a number of 
physical properties of SNe which we can determine from radio observations and modeling. 

2.4.1. Mass-Loss Rate from Radio Absorption 

From the Chevalier (1982a,b) model, the turn-on of the radio emission for RSNe pro- 
vides a measure of the presupernova mass-loss rate to wind velocity ratio (M/wwind)- For 
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the case in which the external absorption is entirely due to a uniform medium, using the 
formulation of Weiler et al. (1986) Eq. 16, we can write 



1.5 



(w^ind/10 km s"^) csM,„.fo,^ VlO^kms"^ 

^ N 1.5 / ^ N 1.5m / rp \ 0-68 



45 days/ \tj \10* K 



(11) 



Since the appearance of optical lines for measuring SN ejecta velocities is often delayed a 
bit relative to the time of the explosion, we arbitrarily take ti — 45 days. Because our 
observations have shown that, generally, 0.8 < m < 1.0 and from Eq. 11 M oc ^ 
the dependence of the calculated mass-loss rate on the date of the initial ejecta velocity 
measurement is weak, M oc tf^'^. Thus, we can take the best optical or VLBI velocity 
measurements available without worrying about the deviation of their exact measurement 
epoch from the assumed 45 days after explosion. For convenience, and because many SN 
measurements indicate velocities of ~ 10, 000 km s~^ , we take Vi = Wbiastwave = 10, 000 
km s~^ . We also normally adopt values of T = 20, 000 K, Wwind ~ 10 km s~^ (which is 
appropriate for a RSG wind) and, generally from our best fits to the radio data for each 
RSN, t = (teem peak " ^o) days, 

'^csMuniform fro^i Eq. 3, aud m from Eq. 6 or 7, as appropriate. 



2.4-2. dumpiness of the Presupernova Wind 

In their study of the radio emission from SN1986J, Weiler, Panagia, & Sramek (1990) 
found that the simple Chevaher (1982a,b) model could not describe the relatively slow turn- 
on. They therefore added terms described mathematically by TcsMdumps ^Iqs. 1 and 5. This 
extension greatly improved the quality of the fit and was interpreted by Weiler, Panagia, 
& Sramek (1990) to represent the possible presence of filamentation or dumpiness in the 
GSM. 

Such a dumpiness in the wind material was again required for modeling the radio data 
from SN1988Z (Van Dyk et al. 1993) and SN1993J (Van Dyk et al. 1994). Since that time, 
evidence for filamentation in the envelopes of SNe has also been found from optical and UV 
observations (e.g., Filippenko, Matheson, & Earth 1994; Spyromilio 1994). 

In the presence of a clumpy absorbing medium, the expression to derive the mass- 
loss rate given by Eq. 11 needs to be revised. For the case of S = S', it is clear that the 
fraction of clumpy material remains constant throughout the whole wind established CSM 
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and, therefore, the radio signal from the SN suffers an absorption rcsMunifo^^ from the uniform 
component of the CSM plus an additional absorption with an even probability distribution 
between and TcsMdumps fro™ the clumpy or filamentary component of the CSM. Since the 
mass-loss rate is proportional to the square root of the optical depth, it will be enough to 
replace r^gj^j in Eq. 11 with the square root of an effective optical depth, t^q (see also 
Van Dyk et al. 1994), which is the appropriate average over the possible extremes of the 
optical depth and is defined as 

< T-eff^ >= 0.67 [(TcSM„„ifo,^ + TCSM,i,^pJ^'^ - '^CSM,„ifo,^]T"cSMeiumps (^2) 

SO that Eq.(ll) can more generally be written as 



(wwind/10 km s"^) ■ VlO^kms 



45 days/ \tj K 



(13) 



which, in the limit of tcsm i ^ — approaches Eq. 11 and in the limit of tcsm „if — 



approaches Eq. 11 with TcIm replaced by 0.67 Tc|m, 



clumps 

uniform ^ -'■^■'-clumps 

While intermediate cases where b ^ S', TcsMu„ifo„n 0) ''"csMdumps will yield results 
with larger errors, it is felt, considering other uncertainties in the assumptions, that Eqs. 11 
and 13 yield reasonable estimates of the mass-loss-rates of the presupernova star. Absorption 
mass-loss rate estimates from Eq. 11 or 13, as appropriate, for Type Ib/c SNe tend to be 
10-6 M© yr-i while those for Type II SNe tend to be higher at ~ 10"^ - 10"^ M© yr'^ 



2.4-3. Mass-Loss Rate from Radio Emission 

For comparison purposes, we can also try to estimate the presupernova mass-loss rate 
which established the CSM by considering the radio emission directly. Prom the Chevaher 
model and Weiler et al. (1989) one can write the radio emission from an RSN in the form: 



M(Mq yr ^) ^ Q ^ 1 Q-9 f -^6cm peak A "'^^ / ^ecm peak - V'"^^ , . 

.ind/lOkms-^) ■ VlO^'ergss-iHz-i; ^ (days) J ^ ^ 



(wwind/lO km s ) 
for Type Ib/c SNe and 
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M(Mq yr ^) ^ 2 X 10"*^ ( -^6cm peak \ "'^^ / ^6cm peak - ^0 \ 

Kind/10 km s-^) ■ V 10^6 ergs s-i Hz- V V (days) J ^ ' 

for Type II SNe, assuming that the absorption Teem peak ~ 1, from whatever origin, at the 
time of observed peak in the 6 cm flux density. 

The coefficients in Eqs. 14 and 15 depend on the amount of kinetic energy that is trans- 
ferred to accelerate relativistic electrons and on the details of the acceleration mechanism. 
Although Fermi acceleration is generally accepted as the relativistic electron acceleration 
process and it is usually assumed that some fixed fraction of the explosion kinetic energy 
is transformed into relativistic synchrotron electrons (often assumed the physics of 

these two aspects is not known in detail a 'priori. Therefore, Eqs. 14 and 15 can only be 
"calibrated" by using the values of well studied RSNe and the assumption that all RSNe of 
the same type will have similar characteristics. The constants in Eqs. 14 and 15 have thus 
been determined from averages within the two RSN subtypes (Type Ib/c and Type II) of 
those RSNe which have pure, uniform absorption {i.e., Kj, — 0), with the omission of SN 
1987A because of its blue supergiant (BSG) rather than red supergiant (RSG) progenitor. 

It should be kept in mind that, because the detailed mechanism of radio emission from 
SNe is not well understood and the estimates have to rely on ad hoc calibrations from the 
few RSNe which have well measured hght curves, mass-loss rate estimated from Eqs. 14 and 
15 can only complement, and perhaps support, the more accurate determinations done from 
radio absorption using Eqs. 11 and 13. 



2.1^.1^. Changes in Mass-loss Rate 

A particularly interesting case of mass-loss from an RSN is SN1993J, where detailed 
radio observations are available starting only a few days after explosion. Van Dyk et al. 
(1994) find evidence for a changing mass-loss rate for the presupernova star which was 
as high as ~ 10"^ M0 yr"^ approximately 1000 years before explosion and decreased to 
~ 10"^ M0 yr"^ just before explosion, resulting in a relatively flat density profile of p oc r~^-^. 

Pransson & Bjorgsson (1998) have suggested that the observed behavior of the f-f 
absorption for SN1993J could alternatively be explained in terms of a systematic decrease 
of the electron temperature in the circumstellar material as the SN expands. It is not 
clear, however, what the physical process is which determines why such a cooling might 
occur efficiently in SN1993J, but not in SNe such as SN1979C and SN1980K where no such 
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behavior is required to explain the observed radio turn-on characteristics. Also, recent X-ray 
observations with ROSAT of SN 1993J imply a non-r~^ CSM density surrounding the SN 
progenitor (Immler, Aschenbach, & Wang 2001), with a density gradient of p oc r~^^^ . 

Moreover, changes in presupernova mass-loss rates are not unusual. Montes et al. 
(2000) find that Type II SN1979C had a slow increase in its radio hght curve after day 
~ 4300 which implied an excess in flux density by a factor of ~1.7 with respect to the stan- 
dard model, or a density enhancement of ~ 30% over the expected density at that radius. On 
the other hand. Type II SN1980K showed a steep decline in flux density at all wavelengths 
by a factor of ~ 2 occurring between day ~3700 and day ~4900. Such a sharp decline in flux 
density implies a decrease in pcsM by a factor of ~ 1.6 below that expected for a r^^ CSM 
density profile (Montes et al. 1998). If one assumes the radio emission arises from a ~ 10^ 
km s~^ blastwave travehng through a CSM established by a ~ 10 km pre-explosion stel- 
lar wind, this implies a significant change in the stellar mass-loss rate, for a constant speed 
wind, ~ 12, 000 yr before explosion for both SNe. 

3. SN1998bw 

The radio emission from SN 1998bw reached an unusually high 6 cm spectral luminosity 
at peak of ~ 6.7 x 10^^ erg s"^ Hz~^, i.e., ~ 3 times higher than either of the well studied, 
very radio luminous Type Iln SNe SN 1986J and SN 1988Z, and ~ 40 times higher than 
the average peak 6 cm spectral luminosity of Type Ib/c SNe. It also reached this 6 cm peak 
rather quickly, only ~ 13 days after explosion. 

SN 1998bw is unusual in its radio emission, but not extreme. For example, the time 
from explosion to peak 6 cm luminosity for both SN 1987A and SN 1983N was shorter and, 
in spite of the fact that SN1998bw has been called "the most luminous radio supernova ever 
observed," its 6 cm spectral luminosity at peak is exceeded by that of SN 1982aa (Yin 1994). 
However, SN 1998bw is certainly the most radio luminous Type Ib/c RSN observed so far 
by a factor of ~ 25 and it reached this higher radio luminosity very early. 

3.1. Expansion Velocity 

Although unique in neither the speed of radio light curve evolution nor in peak 6 cm 
radio luminosity, SN1998bw is certainly unusual in the combination of these two factors - 
very radio luminous very soon after explosion. Kulkarni et al. (1998a,b) have used these 
observed qualities, together with the lack of interstellar scintillation (ISS) at early times. 
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brightness temperature estimates, and physical arguments to conclude that the blastwave 
from SN1998bw giving rise to the radio emission must have been expanding relativistically. 
On the other hand, (Waxman & Loeb 1999) argue that a sub-relativistic blastwave can 
generate the observed radio emission. However, both sets of authors agree that a very 
high expansion velocity (;^ 0.3c) is required for the radio emitting region under a spherical 
geometry. 

Simple arguments confirm this high velocity. To avoid the well known Compton Catas- 
trophe, Kellermann & Pauliny-Toth (1969) have shown that the brightness temperature 
Tb < 10^^ K must hold and Readhead (1994) has better defined this limit to Tb < 10^^'^ 
K. From geometrical arguments, such a limit requires the radiosphere of SN1998bw to have 
expanded at an apparent speed ^ 230, 000 km s~^ , at least during the first few days after 
explosion. While such a value is only mildly relativistic (F ~ 1.6; F = 1^2 )^ ^^^^^ 

unusually high. Measurements by Gaensler et al. (1997) and Manchester (2001) have 
demonstrated that the radio emitting regions of the Type II SN1987A expanded at an aver- 
age velocity of ~ 35, 000 km s~^ over the 3 years from 1987 February to mid- 1990 so that, in 
a very low density environment such as one finds around Type Ib/c SNe, very high blastwave 
velocities appear to be possible. 



3.2. Radio Light Curves 

An obvious comparison of SN1998bw with other RSNe is the evolution of its radio flux 
density at multiple frequencies and its description by known RSN models. The available 
radio data can be found on 

http://www.narrabri. atnf. csiro. au/ mwiermg/grh/grh980425/ 

and are plotted in Fig. 2. SN1998bw shows an early peak which reaches a maximum as early 
as day 10 - 12 at 8.64 GHz, a minimum almost simultaneously for the higher frequencies 
(z/ > 2.5 GHz) at day ~ 20 - 24, then a secondary, somewhat lower peak at later times after 
the first dip. An interesting characteristic of this "double humped" structure is that it dies 
out at lower frequencies and is relatively inconspicuous in the 1.38 GHz radio measurements 
(see Fig. 2). 

Li & Chevalier (1999) propose an initially synchrotron self-absorbed (SSA), rapidly 
expanding blastwave in a p oc r^^ circumstellar wind model to describe the radio light curve 
for SN 1998bw. This is in many ways similar to the Chevaher (1998) model for Type Ib/c 
SNe, which also included SSA. To produce the second peak in the radio light curves Li & 
Chevaher (1999) postulate a boost of blastwave energy by a factor of ~ 2.8 on day ~22 in 
the observer's time frame. 
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Modeling of the radio data for SN 1998bw with the well estabhshed formahsm for RSNe 
presented above shows that such an energy boost is not needed. A fast shock interacting 
with a dense, slow, stellar wind established, ionized CSM, which is modulated in density 
over time scales similar to those seen for RSNe, can produce a superior fit to the data. 
No blastwavc reacceleration is required and no synchrotron self- absorption (SSA) at early 
times is apparent. The parameters of our best fit model are given in Table 1 and shown 
as the curves in Fig. 2. Li & Chevaher (1999) do not give a quantitative indication of the 
quality of their fit. However, a visual comparison of the curves in Fig. 2 with those of Li & 
Chevaher (1999) Fig. 9, shows that our purely thermal absorption model with structured 
CSM provides a superior fit. 



4. Discussion 

4.1. Radio Light Curve Description 

An unusual feature of the SN1998bw radio light curves shown in Fig. 2 is its apparently 
"double humped" structure with a second peak in the radio emission near Day ~ 30. Such 
a "double humped" structure of the radio light curves can be reproduced by a single energy 
blastwave encountering differing CSM density regimes as it travels rapidly outward. This is a 
reasonable assumption since complex density structure in the CSM surrounding SNe, giving 
rise to structure in the radio light curves, is very well known in such objects as SN1979C 
(Weiler et al. 1991, 1992a; Montes et al. 2000), SN1980K (Weiler et al. 1992b; Montes et 
al. 1998), SN1988Z (Van Dyk et al. 1993; Lacey et al. 2001) and, particularly, SN1987A 
(Jakobsen et al. 1991). 

Additionally, what has not been previously recognized is a sharp drop in the radio 
emission near Day ~ 75 and a single measurement epoch at Day 192 which is significantly 
(~ 60%) higher at all frequencies than expected from the preceeding data on Day 149 and 
the following data on Day 249. 

Quite surprisingly, both of these temporary increases in radio emission can be explained 
by the SN blastwave encountering physically similar shells of enhanced density. The first 
enhancement or "bump" after the initial outburst peak is estimated to start on Day 25 and 
end on Day 75, i.e., having a duration of ~ 50 days and turn-on and turn-off times of about 
12 days, where the increased radio emission and absorption {Ki and K^) by a factor of ~ 1.6 
imply a density enhancement of ~ 30%. Exactly the same density enhancement factor and 
length of enhancement is compatible the "bump" observed in the radio emission at Day 192 
{i.e., the single measurement within the 100 day gap between measurements on Day 149 and 
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Day 249), even though the logarithmic time scale of Fig. 2 makes the time interval look much 
shorter. Unfortunately, the decreased sampling interval has only one set of measurements 
altered by the proposed Day 192 enhancement, so that we cannot determine its length other 
than being < 100 days. 

Note that the fit hsted in Table 1 and shown as the curves in Fig. 2 requires no "uniform" 
absorption {K2 = 0) so that all of the f-f absorption is due to a clumpy medium as described 
in Eqs. 1 and 5. These results, combined with the estimate of a high blastwave velocity, 
suggest that the CSM around SN1998bw is highly structured with little, if any, inter-clump 
gas. The clump filling factor has to be high enough to intercept a considerable fraction of the 
blastwave energy and low enough to let radiation escape from any given clump without being 
appreciably absorbed by any other clump, i.e., a filling factor of ~0.1 - 0.3. The blastwave 
can then easily move at a speed which is a significant fraction of the speed of light, because 
it is moving in a very low density medium, but still cause strong energy dissipation and 
relativistic electron acceleration at the clump surfaces facing the SN explosion center. 

If the knots are dense and relatively opaque to radio emission, we mainly observe radia- 
tion produced at the surfaces of the knots in the CSM on the far-side of the SN in which some 
internal absorption may also be occurring. The radiation from near-side knots is probably 
absorbed internally by the knots themselves and, therefore, lost from the signal we detect. 

Note also from the fit given in Table 1 that the presence of a K4 factor implies there 
is a more distant, uniform screen of ionized gas surrounding the exploding system which is 
too far to be affected by the rapidly expanding blastwave and provides a time independent 
absorption. 

4.2. Physical Parameter Estimates 

Using the fitting parameters from Table 1 and Eq. 13, we can estimate a mass-loss rate 
from the pre-explosion star. However, the proper parameter assumptions are rather uncertain 
for these enigmatic objects. For a preliminary estimate we shall assume ti = 23 days, 
t = (^ecm peak " ^o) = 13.3 days, m = -{a - /? - 3)/3 = 0.78 (Eq. 7), unwind = 10 km s"^ (for 
an assumed RSG progenitor), Vi = Vbiastwavc = 230, 000 km s^^ , and T = 20, 000 K. We shall 
also assume, since the radio emission implies that the CSM is highly clumped {i.e., K2 = 0), 
a filling factor, /, of only / ~ 10%, which reduces the estimated mass loss rate by a factor of 
y/f. Within these rather uncertain assumptions, Eq. 13 yields an estimated mass-loss rate 
of M ~ 3.5 X 10~^ M0 yr~^ with density enhancements of ~ -\/L6 or ~ 30% during the two 
known, extended "bump" periods. 
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Assuming that the blastwave is travehng at a constant speed of ~ 230, 000 km s'^ the 

start of the first "bump" on Day 25 imphcs that it starts at ~ 5.0 x 10^^ cm and ends on 
Day 75 at ~ 1.5 x 10^^ cm from the star. Correspondingly, if it was estabhshed by a 10 
km s^^ RSG wind, the 50 days of enhanced mass-loss ended ~ 1, 600 yr and started ~ 4, 700 
yr before the explosion. The earlier high mass-loss rate epoch indicated by the enhanced 
emission on Day 192 in the measurement gap between Day 149 and Day 249 implies, with 
the same assumptions, that it occurred in the interval between ~ 9, 400 yr and ~ 15, 700 yr 
before explosion. It is interesting to note that the time between the presumed centers of the 
first and second increased mass-loss episodes of ~ 9,400 yr is comparable to the ~ 12,000 
yr before explosion at which SN1979C had a significant mass-loss rate increase (Montes et 
al. 2000) and SN1980K had a significant mass-loss rate decrease (Montes et al. 1998), 
thus establishing a possible characteristic time scale of ~ lO"' yr for significant changes in 
mass-loss rate for pre-explosion massive stars. As discussed by Panagia & Bono (2001) 
and Bono & Panagia (2001), mass-loss variations with such time scales can be explained in 
terms of pulsational instabilities of red supergiants with initial masses of ~ 20 M©. 



5. GRB Radio Afterglows and Type Ib/c RSNe 

Since the suggestion of a possible relation between SN1998bw and GRB980425, it has 
remained a tantalizing possibility that the origin of at least some GRBs is in the better 
known Type Ib/c SN phenomenon. First, of course, one must keep in mind that there may 
be (and probably is) more than one origin for GRBs, a situation which is true for most other 
classes of objects. For example, SNe, after having been identified as a new phenomenon in 
the early part of the last century, quickly split into several subgroups such as Zwicky's Types 
I - VI, then coalesced back into just two subgroupings based on Ha absent (Type I) or Ha 
present (Type II) in their optical spectra. This simplification has not withstood the test of 
time, however, and subgroupings of Type la, lb, Ic, II, lib, Ilpec, Iln, and others have come 
into use over the past 20 years. 

GRBs, although at a much earfier stage of understanding, have similarly started to split 
into subgroupings. The two currently accepted groupings are referred to as "fast-hard" and 
"slow-soft" from the tendency of the 7-ray emission for some to evolve more rapidly (mean 
duration ~0.2 s) and to have a somewhat harder spectrum than for others which evolve more 
slowly (mean duration ~20 s) with a somewhat softer spectrum (Fishman & Meegan 1995). 

Since we are only concerned with the radio afterglows of GRBs here, all of our examples 
fall into the "slow-soft" classification, at least partly because the "fast-hard" GRBs fade 
too quickly for foUowup observations to obtain the precise positional information needed 
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for identification at longer wavelengths. It is therefore uncertain at present whether "fast- 
hard" GRBs have radio afterglows or even whether the "slow-soft" GRBs represent a single 
phenomenon. If, however, we assume that all "slow-soft" GRBs have a similar origin and 
that GRB980425/SN1998bw is a key to this puzzle telhng us that "slow-soft" GRBs have 
their origin in Type Ib/c SNe, we can investigate relations between the two observational 
phenomena. 

Keeping in mind that much of the radio data available on Type Ib/c SNe is of poor 
quality due to their relative faintness, and that the situation is even worse for GRB radio 
afterglows, some of the similarities are: 

1. GRB980425 and Type lb SN1998bw are possibly related. 

2. Both GRBs and Type Ib/c SNe are suggested by some models to originate from the 
core collapse of a massive (perhaps stripped) progenitor star. 

3. Both GRBs and Type Ib/c SNe are suggested, by some models, to release ~ 10^^ — 10^^ 
ergs of kinetic energy in a few seconds. 

4. Both GRBs and Type Ib/c SNe are detectable at cm wavelengths very soon after the 
explosion - within a few days, while Type II SNe generally have much longer delays. 

5. Both GRBs and Type Ib/c RSNe generally have relatively steep radio spectra (a ^ 
— 1; S* oc z/"*"") compared to, for example, extragalactic QSOs or radio galaxy nuclei 
{a ;^ -0.75) or Type II RSNe {a ;^ -1). 

There is presently one clear difference between known GRBs and Type Ib/c SNe, 
however. Type Ib/c RSNe can generally be treated as classical, non-relativistic phenom- 
ena, sufficiently close that cosmological effects can be neglected, while most known GRBs 
appear to be highly relativistic phenomena (F ~ 10) at cosmological distances (z ~ 1). 
SN1998bw/GRB980425, which is not very distant (z = 0.0085, ~40 Mpc) and only mildly 
relativistic (F ~ 1.6), is possibly the "bridge" between the two groups. We shall discuss the 
properties of specific GRBs in more detail in a following paper. 

While observational similarity does not necessarily mean identical physical processes, 
the possibihty of similar or related phenomena giving rise to both GRBs and Type Ib/c 
SNe must be considered. Speculatively, GRBs may represent a rare subclass of Type Ib/c 
SNe and/or the subset of Type Ib/c SNe where we he within the narrow emission cone of 
a relativistic jet. In either case, they are sufficiently rare that they can only be found by 
sampling very large volumes of space, i.e., preferentially at high redshifts. 
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6. Summciry 

We have collected and rationalized a parameterization and physical model description 
which has been very successful in explaining the radio emission from SNe. This we have then 
applied to SN1998bw, which is thought to be related to GRB980425. From this parameter- 
ization of the radio emission for SN1998bw, we are able to describe the radio light curves 
for SN1998bw in previously unattainable detail and show that a mildly relativistic blastwave 
interacting with a clumpy, structured CSM established by the wind from a massive progen- 
itor provides a satisfactory description. Superimposed on a relatively constant mass-loss of 
~ 3.5 X 10~^ M0 yr~^ are epochs of higher mass-loss rate which increase the CSM density by 
a factor of ~ 30%. At least two of these episodes can be identified and occur at an interval 
of ~ 9, 400 yr and a duration of ~ 3, 000 yr. 

KWW & MJM wish to thank the Office of Naval Research (ONR) for the 6.1 funding 
supporting this research. Additional information and data on radio emission from SNe and 
GRBs can be found on http://rsd-www.nrl.navy. 'mil/72 14/weiler/ and linked pages. 
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Table 1. SN1998bw Modeling Results. 



Parameter 



Value 



K2 
S 

S' 



a 



-0.71 
-1.38 
2.4 X 10^ 




1.7 X 10^ 

-2.80 
1.2 X 10-2 



to (Explosion Date) 1998 Apr. 25.90915 



^Enhanced by a factor of 1.6 over the intervals 
Day 25 - 75 and Day 165 - 215, although the latter 
interval could be as long as 100 days and still be 
compatible with the available data. 

Assuming ti = 23 days, t = (teem peak - to) = 
13.3 days, m = —{a — P — 3)/3 — 0.78, Wwmd — 
10 km , Vi = ^^biastwave = 230,000 km s"^ , 
T = 20, 000 K, and filling factor, / = 0.1. 




13.3 

37.4 

38.9 
6.7 X 1028 
3.5 X 10-5 
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Fig. 1. — Cartoon, not to scale, of the SN and its shocks along with the stellar wind 
established CSM and more distant ionized material. The radio emission is thought to arise 
near the blastwave. The expected locations of the several absorbing terms in Eqs. 1-10 are 
illustrated. 

Fig. 2. — The radio hght curves of SN 1998bw at 8.64 GHz (3.5 cm; upper left, open circles, 
solid line), 4.80 GHz (6.3 cm; upper right, stars, dashed line), 2.50 GHz (12 cm; lower left, 
open squares, dash-dot line) and 1.38 GHz (21.7 cm; lower right, open triangles, dotted line). 
The curves are derived from a best fit model described by the Eqs. 1-10 and the parameters 
and assumptions listed in Table 1. During the 50 day intervals from Day 25 - 75 and from 
Day 165 - 215 the emission and absorption factors {Ki and K^) are increased by a factor of 
1.6 with a 6 day boxcar smoothed turn-on and turn-off of the enhanced emission/absorption. 
The physical interpretation of these enhancements is discussed in §4. 



